Electronic fiber in graphene 
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We investigate theoretically the transmission properties through a p-n-p junction on graphene. 
Here we show that the electronic transport property presents deep analogies with light propagation. 
It originates from the similarity between the linear spectrum of the Dirac fermions and photons 
that obey the Maxwell's equations. We demonstrate the p-n-p channel acts as an electronic fiber in 
which electrons propagate along the channel without dissipation. 
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In recent years, graphene, a single layer of carbon atoms 
arranged in a hexagonal lattice, exhibits abundant physics 
and potential applications. [1-5] Quantum transport prop- 
erties in graphene have attracted increasing attentions both 
from the fundamental physics and potential application in 
carbon-based nano-electronics devices. The properties arise 
from the unique linear energy dispersion and the chiral nature 
of electrons at the K (K ) point of the Brillouin zone, the en- 
ergy spectrum exhibits a linear dispersion that can be well 
described by the massless Dirac equation. [5[ These massless 
Dirac fermions, can be viewed as electrons that lose their rest 
mass mo, or as neutrinos that acquire the electron charge e. 
Electrons in graphene are quite different from the electrons in 
conventional semiconductor two dimensional electron gases. 

Dirac equation describing the motion of massless quasi- 
particles is mathematically similar to the Helmholz equation 
for an electromagnetic wave. Graphene nano structures have 
been achieved by lithography as small as few nanometers. [6|] 
While, the mean- free path of electron in graphene approaches 
the order of micron at room temperature, [4] and the electron 
wavelength is even lager. [7( In such ballistic regime, the scat- 
tering of electrons by potential barriers can be understood by 
comparing with the reflection, refraction and transmission of 
electromagnetic waves in inhomogeneous media. So, graphene 
monolayer is a suitable candidate to examine the optical-like 
phenomena of the Dirac fermions, that have attracted a lot of 
interest, gun such as waveguide, Goos-Hanchen effect, [13| ] 
Bragg reflection, coherent buffers/memories, and Brewster- 
type angle in graphene. Note that it was clearly demon- 
strated by Cheianov et al. [8] that graphene based n-p-n junc- 
tions could be used as electronic lenses or beam splitters. In 
this letter, we investigated theoretically transport property of 
Dirac electrons in a p-n-p channel created by two electrostatic 
gates that are deposited on top of the graphene monolayer (see 
Fig. [Tj). The transport process along the channel determined 
by the multiple total internal reflection. We find that in such 
a p-n-p channel, electrons can be confined in between two in- 
terfaces, in analogy with the light propagation in an optical 
fiber. In addition, the approach presented in this work could 
also be applied directly to electronic excitations in another 
recently hotspot material, namely the topological insulator 
that also has Dirac cones. [l7L [l8j] 

Model — Tunable potential barriers could be created on a 
single- layer graphene sheet, [l9L [20(] by electrostatic gates de- 



posited on top or back of the graphene monolayer. The elec- 
trostatic potentials can be tuned by changing applied voltages 
(see Fig.[T]). The low-energy electrons near the K point of the 
Dirac cones can be well described by the effective Hamilto- 
nian [2ll] 

H = hv F cr-k + V i (r,t), (1) 

where the superscript i indicates the i-th region, vf is the 
Fermi velocity, a is the Pauli matrix and V % is the height 
of the electrostatic barrier. The wavevector k satisfies k t2 + 
k y 2 = (E — V 1 ) 2 1 \Tivf) 2 '. Note that the translational invari- 
ance along the y direction gives rise to the conservation of k y , 
and thus the solutions can be written as i/)(x, y) — ip{x)e lkyV . 
For brevity, we will set Tivf = 1 in the follows. 

Refraction and total reflection — We started by investi- 
gating electron tunneling through an electric p-n junction 
on graphene, with V tn = in the n-doped region and 
yout _ y j n p_d p ec l region. This requirement fixes 
the dependence of t he longitudinal out going wave vector 
k' x = sgn(E F - V)y / (E F - V) 2 - (k y ) 2 . Note that, when 
Ef — V < 0, sgn(EF — V) = —1 denotes the interband 
scattering at the interface of the p-n junction which is re- 
quired by the helicity nature of the Dirac Fermions. Conti- 
nuity of the wave functions at x = gives the transmission 
coefficient: t = 2k x (E F - V)/(E F k x + E F k x - k x V + ik y V). 
Since the modes in n-doped and p-doped regions have differ- 
ent group velocities, the transmission probability is given by 

|2 „,< _ 1 djEF-V*) _ K 
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FIG. 1: Schematic of a p-n-p channel in graphene. The shaded re- 
gions denote the regions below the electric gates. The lower pannel 
describes the potential profile. 
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FIG. 2: (a) Transmission probability as a function of the incident 
angle for several representative gate voltage V ranging from 30 meV 
to 120 meV. (b) Transmission probability as a function of the gate 
voltage for different incident angles. The incident energy is fixed 
at Ep = 80 meV both in panel (a) and (b). (c) The phase diagram 
in momentum space. The outer circle indicates the wavevector 
in the incident n-doped region and the inner circle indicates the 
wavevector in the transmitted p-doped region. 

gating states and v % — for evanescent states. This expres- 
sion has the advantage that can well describe the propagating 
states and evanescent states in a unified form. We can also 
define the incident and refractive angles for electrons in propa- 
gating states (fin = arctan(/c y //c x ), and ip ou t = arctan(/c^//c x ), 
and rewrite the transmission probability as 

T= 4cos((^ n ) 2 

v in ' 2 + 2sgn(E F - V) cos(<^w + <p in + S[-(E F - V)} ■ tt) ' 

(2) 

This is valid for incident angle Lpi n < ip c = arcsin(|£ , F — 
V\/Ef), and the Heaviside step function G guarantees a gen- 
eral expression of the transmission probability whether in- 
terband scattering occurs or not. ip c is analogy to critical 
angle for total reflection in optics and sm(ip ou t)/ sm(<^ n ) = 
Ef/(Ef — V) = n gives the Snell's law for transmitted elec- 
trons. It shows that the effective Fermi energy (Ef — V) plays 
the role of index of refraction in optical medium. Note that 
when Ef — V < 0, refractive index n of graphene is negative 
as a metamaterial. In this case, electrons are injected to the 
hole branch in p region, the chirality of Dirac electrons in 
graphene gives rise to sign reversal of the momentum of the 
forward-going state. Importantly, one can tune the refractive 
index n in quite a large range via electrostatic gates, which is 
not easy for normal metamaterials. We note en passant that 
such optic-like behaviours are absent in the tunneling through 
a barrier, [22] since the effective Fermi energies at both side of 
the barrier are Ef, indicating an identical index of refraction. 

The transmission spectra for electrons traversing such a 
p-n junction are shown in Fig. [2ja) and (b), which can be 
understood from the phase diagram of Fig. EJc). For a fixed 
incident energy, the transmission declines sharply and then is 
blocked when the incident angle ipi n exceeds a critical value 
Lp c (see Fig. [2ja)), while the wavevector k x in the p-doped 
region become imaginary denoting the appearance of evanes- 
cent modes. The critical angle cp c for total reflection can be 
estimated by the Snell's law, cp c = arcsind^F — V\/Ef)- The 
transmission can also be tuned by the gate voltage, or equiva- 
lently, the barrier height V (see Fig.|2jb)). A transmission gap 
appears as the gate voltage V approaching the incident energy 
Ef and become wider as increasing the incident angle ifi n . 
These tunneling features are well appreciated by inspecting 
the phase diagram as shown in Fig.[2jc). The Fermi wavevec- 



tor is different in n-doped and p-doped region as a conse- 
quence of gate voltage V, i.e., kF,N = Ef, kF,p — \Ef — V\, 
while the component k y = kF,N sin (ipi n ) is conserved. The 
requirement k y < kF,p leads to non-zero transmission for al- 
lowed incident angles as a sector in Fig. [2jc) in accordance 
with our numerical results in Fig. [2ja). It is clear that kF,p 
declines as the gate voltage V approaching the incident en- 
ergy Ef, so transmission gaps appear when kF,p drops below 
k y that verifies the calculation as shown in Fig. [2jb). 
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FIG. 3: Energy spectrum for the few lowest channel modes for 
V = 160raeV, D = 80nra. The inset shows the probability for 
states indicated by solid circular dot with the same color. 
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FIG. 4: (a) Conductance versus incident energy. TV labels the 
number of atoms contained in a cell of the middle channel, (b)-(h) 
Density distribution of electron states corresponding solid circular 
dot in Panel (a) for (b)-(f) N = 48, V = 0.8 eV, and (g)-(h) 
N = 20,V = 1.6 eV. 

Fiber guiding-We considered an electric p-n-p channel with 
width D (see FigHJ and focus on the electrons transporting 
along the interface rather than across it. The transmission 
gap discussed above can lead to confinement of electrons be- 
tween the two interfaces associated with multiple total inter- 
nal reflections. Electrons can transport along the p-n-p chan- 
nel just like photons in an optical fiber. We calculated the 
energy dispersion relation of the bound states in the case of 
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total internal reflection at both two interfaces (i.e. ipi n > ip c , 
k' x is imaginary). Matching of propagating waves to evanes- 
cent waves at x = — D and x = 0, the energy spectrum of 
the bound states can be obtained from the zero point of the 
determinant of the coefficients, which can be reduced to a 
transcendental equation — ik' x k x cos(k x D) + [ky — Ef(Ef — 
V)] s'm(k x D) = 0. 

The energy spectrum Ep{k y ) is plotted for the first few 
bound states in Fig. [3] The dashed lines \k y \ — \Ef — V\, and 
\k y \ = \Ef\ give the boundary that delimit the continuum re- 
gion, corresponding to propagating states in ^-direction. Con- 
finement of electrons in ^-direction gives rise to non-dispersive 
transporting in ^/-direction. The bound states between the 
two interfaces thus play the roles of guide modes along the 
channel that can effectively convey the charge carriers with 
group velocities indicated by the slope of the dispersion rela- 
tion. The spectrum is symmetric with respect to the trans- 
verse wave vector k y = 0, guaranteed by the time reversal 
symmetry. So the transport properties is isotropic for each 
incident terminal. The inset of Fig. [3] shows the probability 
density of the bound states marked by the solid circular dot. 
The vertical dashed lines indicate the interfaces of the p-n-p 
junction. It is clear that the solutions are standing waves in 
the channel and decay exponentially in the barrier regions. 
Note that the order of the bound states equals to the number 
of peaks of standing waves in the channel, e.g. the 1st, 3rd 
and 5th bound states or guide modes are shown in the inset 
of Fig. [3] Thus, the bound states serve as guide modes in 
cavity. Finally, we focused on how the fiber guiding features 
shown above are reflected in the measurable quantity, the con- 
ductance G. We performed numerical simulations of electri- 
cal conduction in a tight-binding model of a graphene zigzag 
ribbon covered by electrodes. The conductance is evaluated 
using the recursive Green function technique and Landauer- 
Biittiker formula. The conductance of the gated ribbon dis- 
plays a step-like feature, which corresponds to the opening of 
the new modes as the Fermi energy increases (see Fig. H{a)). 
When Ef is very small, there is no bound states in the chan- 



nel, the conductance is contributed by the edge states of the 
zigzag ribbon, since our simulation are performed in a lattice 
scale. In such low incident energy region, the conductance 
for a thinner channel may also depend on the former hole- 
branch bulk states that are lifted up to the Fermi energy by 
the gate voltage. A lager gate voltage can make this effect 
more pronounced. When the incident Fermi energy lays in 
an interval around the gate voltage, electrons could be to- 
tally reflected with a small enough critical angle since there 
exist a great difference between the refraction indexes at two 
sides of the interface. In this case, most transmission modes 
are formed by the bound states in the channel and desired 
high effective fiber guiding is obtained. As the Fermi energy 
is increased continuously, the refraction indexes of the chan- 
nel and the gated region approach to each other gradually, so 
the total reflection critical angle becomes very large. Elec- 
trons can leak out of the channel easily, leading to decrease of 
guide efficiency. This is because the transmission modes are 
mixture of both bound states and bulk states out of the chan- 
nel. Above expectations are demonstrated by examining the 
electron density distributions as shown in Fig. [UJb)-(h)). It is 
apparent that electrons are mainly distributed in the channel 
when the incident Fermi energy is around the gate voltage and 
thus the conductance depends on the guide modes by virtue 
of the gate voltage and channel width. In addition, similar 
results could be obtained when the gate voltage or incident 
energy are minus due to the particle-hole symmetry. 

In summary, We investigate theoretically the electronic 
quantum transport properties in a p-n-p junction on 
graphene. Our results demonstrate that the proposed struc- 
ture on graphene presents deep analogies with optic phenom- 
ena as refraction and total internal reflection. The quantized 
dispersion relation of the confined electron states are obtained 
in our calculation. We show that these bound states serve as 
guide modes and can carry current along the channel with 
out radiation. The investigation could be helpful to offer a 
functionality of graphene based electronic fiber. 
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